arXiv:1505.07942vl [cond-mat.mtrl-sci] 29 May 2015 


Impact of oxygen doping and oxidation state of iron on the electronic and magnetic 

properties of BaFeOs-^ 

I. V. Maznichenko/’[^ S. Ostanin,^ L. V. Bekenov,^ V. N. Antonov,^’^ 1. Mertig/’^ and A. Ernst^ 

^Institut fiir Physik, Martin-Luther-Universitdt Halle-Wittenberg, D-06099 Halle, Germany 
^Max-Planck-Institut fiir Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany 
^Institute of Metal Physics, Vernadsky Street 36, 03142 Kiev, Ukraine 
(Dated: June 1, 2015) 

We studied structural, electronic and magnetic properties of a cubic perovskite BaFeOa-s (0 < 

5 < 0.5) within the density functional theory using a generalized gradient approximation and a 
GGA+U method. According to our calculations, BaFeOa in its stoichiometric cubic structure should 
be half-metallic and strongly ferromagnetic, with extremely high Curie temperature (Tc) of 700 - 
900 K. However, a such estimate of Tc disagrees with all available experiments, which report that 
Tc of the BaFeOa and undoped BaFeOa-a films varies between 111 K and 235 K or, alternatively, 
that no ferromagnetic order was detected there. Fitting the calculated x-ray magnetic circular 
dichroism spectra to the experimental features seen for BaFeOa, we concluded that the presence of 
oxygen vacancies in our model enables a good agreement. Thus, the relatively low Tc measured 
in BaFeOa can be explained by oxygen vacancies intrinsically presented in the material. Since iron 
species near the O vacancy change their oxidation state from 4-|- to 3-I-, the interaction between 
Fe'*'*' and Fe®'*', which is antiferromagnetic, weakens the effective magnetic interaction in the system, 
which is predominantly ferromagnetic. With increasing 5 in BaFeOa-i, its Tc decreases down to the 
critical value when the magnetic order becomes antiferromagnetic. Our calculations of the electronic 
structure of BaFeOa-a illustrate how the ferromagnetism originates and also how one can keep this 
cubic perovskite robustly ferromagnetic far above the room temperature. 

PACS numbers: 78.20.Ls, 75.50.Bb, 75.30.Mb 


I. INTRODUCTION 

A multiferroic tunnel junction, in which a robustly fer¬ 
roelectric barrier, such as PbTiOa or BaTiOa, is grown 
epitaxially between two ferromagnetic electrodes, allows 
to tune the spin transport by switching the magnetiza¬ 
tion of an electrode. Additionally, the barrier polariza¬ 
tion reversal can change the interfacial magnetoelectric 
coupling and, thus, the four-state tunneling resistance 
can be detected.^ As for the upper electrode, the Fe 
and Co films can be easily grown on the Ti02 terminated 
polar perovskites,!^ while the best bottom-electrode ma¬ 
terial, which can be used here as the substrate, must 
satisfy the following conditions: (i) its structure should 
be in a good lattice match with ferroelectrics and (ii) the 
electrode material should be robustly ferromagnetic with 
relatively high Curie temperature (Tc)- A very good 
choice would be a magnetic perovskite since a ferroelec¬ 
tric can be well grown on such substrates. In particu¬ 
lar, perovskites ABO 3 with a magnetic cation on site B 
are well-known to be compatible with ferroelectrics.^^ Re¬ 
cently, Pantel et alW have demonstrated the four-state 
tunneling resistance effect using the 3.2-nm-thick per¬ 
ovskite barrier PbZro.2Tio,803 (001) grown on ferromag¬ 
netic La(Sr)Mn03 (LSMO). 

One of the possible candidates for the substrate in a 
multiferroic tunnel junction seems to be BaFe03 (BFO). 
Although bulk BFO is hexagonal, the perovskite phase 
can be stabilized in epitaxial films. X-ray diffraction 
studies show that thin BFO films grown on SrTiO3(001) 
adopt the cubic perovskite structure with the lattice pa¬ 


rameter of 3.97 For BaFe03_5 films, the larger lat¬ 
tice parameter of 4.07 A has been measured.® Regarding 
the observed magnetic properties of the cubic BFO, such 
as the Fe magnetic moment, magnetic ordering and Tc, 
the reports seem contradictory. Callender et aZ.® reported 
that a pseudocubic BFO is robustly ferromagnetic (FM) 
with Tc=235 K. Instead, the oxygen deficient BaFe03_5 
films exhibit antiferromagnetism (AFM) at room tem¬ 
perature, while a weakly ferromagnetic behavior extends 
above 390 K.® For thick BFO films, Chakraverty et aZ.® 
reported Tc=115 K and a saturation magnetization of 
3.2 /rs/f.u., which both remain stable as a function of 
film thickness, with no signature of the spin spiral struc¬ 
ture up to 300 nm. 

In the case of completely oxidized stoichiometric cu¬ 
bic BFO, the oxidation state of iron is 4+ that leads to 
the twofold orbitally degenerated configuration t^gCg. In 
polyatomic molecules the half occmied Cg level have to 
split according to Jahn and Teller.l^ In solids the effect 
must be cooperative. However, for AFe03 (A = Ca, Sr, 
Ba) this behavior, known as the cooperative Jahn-Teller 
distortion, was not reported so far. Instead, it occurs 
in antiferromagnetic LaMn03 where up to 750 K the Cg 
splitting opens the band gap.^^i®^ On the other hand, 
the Fe oxidation state in oxides is so unstable that the 
presence of the degenerated Cg level is not required. For 
instance, in metallic CaFe03 the charge disproportion- 
aly is observed: 2Fe^“'' —>■ -|- where 

0 < r < l.ES Thus, the issue of mixed valency might be 
very important in the case of BFO. 

Several theoretical studies of magnetic BFO have been 
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recently made from first principles. For the FM order 
and experimental lattice parameter of the cubic BFO, the 
calculated ma gnetic moment varies between 3.0 hb and 
3.4 per f.u. P^ I ^^ I depending on the model approxima¬ 
tion used to treat the exchange-correlation effects. Ac¬ 
cording to the ab initio investigation by Ribeiro et 
the FM solution in the stoichiometric BFO is energeti¬ 
cally favorable against the possible AFM configurations 
while the FM moment is about 3.6 /is- A helical mag¬ 
netic ordering in BFO, which is usually discussed in the 
context of SrFeOa and CaFe 03 ,^^ has b een investigated 
from first principles by Zhi Li et al.E^llSI xhe authors an¬ 
ticipate that the G-type helical solution can be fixed, 
keeping the reduced volume and, simultaneously, using 
the Hubbard parametrization which treats electronic cor¬ 
relations beyond the standard of the local density ap¬ 
proximation. The theoretical estimations of Tq were not 
performed so far. 

As it was already mentioned above, the Fe atoms in 
the BFO can have various oxidation states. The origin 
of this distinction in real materials can be related for in¬ 
stance with a particular sample preparation procedure, 
that can explain the large discrepancy in results of vari¬ 
ous experiments. In this work, we investigate from first- 
principles the impact of Fe oxidation states on structural, 
magnetic, optical, and electronic properties of chemically 
perfect BFO and BaFeOa-^ (S > 0) with vacancies. The 
paper is organized as follows. In Section we provide 
the essential details of our methods and calculations. In 
Section [ml we present the results and discuss the struc¬ 
tural, electronic, optical, and magnetic properties of the 
BFO. Finally, our summary and conclusions are drawn. 


II. DETAILS OF CALCULATIONS 

Our calculations of the BFO, undoped BaFeOa-a, and 
nominally overdoped BFO were performed using a first- 
principles Green function methocP^I based the multiple¬ 
scattering within the density functional theory in a gen¬ 
eralized gradient approximation (GGA) to the exchange- 
correlation potentially We used a full charge density ap¬ 
proximation to take into account possible non-sphericity 
of the crystal potential and the charge density. 

The Fe oxidation state was modeled using a self¬ 
interaction correction (SIG) methocP^ as it is imple¬ 
mented within the multiple-scattering theory.l^ In this 
case, an oxidation state can be defined by the number of 
valence electrons as 


J^val = Z — {Ncore + Nsic), ( 1 ) 

where Z is the total number of electrons of an isolated 
atom, Nf-ore is the number of electrons occupying the 
deep energy levels and treated as the core electrons, Nsic 
is the number of SIC corrected 3d electrons. 

Since the SIC approach is designed only for ground 
state properties,^ we used a GGA-I-U method^ to take 
into account electronic correlation effects in calculations 



FIG. 1: (Color) The crystalline structure of the cubic BaFeOs. 


of exchange parameters. The later were estimated within 
the magnetic force theorem.^ The critical temperature 
and the magnetic ordering were obtained within the 
Heisenberg model 


H — ^ ^ Jij 7 (2) 

hi 

where i and j label the magnetic atoms, and is a 
unit vector in the direction of the magnetic moment of 
the zth atom. The critical temperatures were estimated 
using both the mean field approximation (MFA) and 
random phase approximation (RPA).^ Both approaches 
were successfully applied to many m agneti cally ordered 
oxides including magnetic perovskites.^yy 

Disorder effects in our study were simulated within 
the coherent-potential approximation (CPA).iyty The 
intrinsic oxygen vacancies in BaFeOs-^ were modeled 
within the same CPA scheme, which implies randomly 
distributed vacancies over all oxygen sites. Regarding 
the Ba and Fe vacancies, in the literature, such cation- 
deficient perovskites are formally denoted as BaFeOa+i, 
although the bulk perovskite structures are unable to ac¬ 
commodate interstitial oxygen. 

The calculations of the XAS and XMCD were per¬ 
formed using the spin-polarized fu lly re lativistic linear- 
muffin-tin-orbital (LMTO) methocpys for the experi¬ 
mentally observed lattice constants. The details are de¬ 
scribed in our previous papers.l^fy 


III. RESULTS 

A. Crystalline structure of BaFeOs 

The ground-state phase of BFO is hexagonal with the 
space group symmetry PV2i/c\^^li contains six formula 
units (f.u.) and two types of environmentally different 
Fe. We calculated the structural properties using both 
the LDAiy and GGA^y to the exchange-correlation po¬ 
tential. The GGA (LDA) yields the equilibrium volume 
of 378.66 A3 (352.84 A^), which is in a good agreement 
with the experimental value of 388.42 A^. For this vol¬ 
ume and chemically perfect BFO, we found that the AFM 
configuration is marginally preferable by 9 meV/f.u., as 
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compared to the FM solution, while the magnetic mo¬ 
ments of each type of Fe differ significantly: 3.34 
and 2.38 fis- The use of the electronic-correlation pa¬ 
rameter U=2 eV applied on iron species increases the Fe 
magnetic moments by 12%, whereas the energy difference 
AE = Efm — Eafm increases up to 28 meV/f.u. 

Concerning the cubic BFO phase, the LDA equilib¬ 
rium lattice parameter of 3.86 A is significantly underes¬ 
timated, as compared to the experimental value of 3.97 A, 
while the use of GGA results in a = 3.963 A. For the lat¬ 
ter we found that the FM ordering is strongly favorable 
by 110 meV/f.u. with respect to the G-type AFM config¬ 
uration. Thus, we suggest that the BFO phase transfor¬ 
mation from its hexagonal AFM to cubic FM structure 
goes with the volume reduction of 1.1% per f.u. and 
the corresponding energy barrier is 0.28 eV/f.u. In BFO 
films, the cubic FM phase can be stabilized at room tem¬ 
perature by anharmonic contributions to the free energy 
and/or the effects of epitaxial strain. 

Thus, the crystalline structure of the hexagonal and 
cubic BFO can be adequately described using a conven¬ 
tional GGA functional. Further, we focus our study on 
electronic, optic, and magnetic properties of the cubic 
BFO phase (see Fig{^. 

B. Electronic and magnetic structure of BaFeOa 

1. GGA approach 

Since we found the GGA to be an appropriate func¬ 
tional in describing the crystalline structure of the cubic 
BFO, we applied the method for further study of its elec¬ 
tronic and magnetic properties. 

Within the GGA, the spin-polarized density of states 
(DOS) of BFO exhibits a metallic behavior in the two 
spin channels, as it is shown in Figj^ This agrees with 
the experimental findings of Gallender et al.^ who ob¬ 
served the ferromagnetic ordering and conductance in 
BFO films, annealed in the oxygen atmosphere. The iron 
atoms appear to have the oxidation state of 2+. There¬ 
with, t 2 g and Cg states are almost occupied in the major¬ 
ity spin channel (4.5 electrons) and only partially occu¬ 
pied (1.5 electrons) in the minority spin channel. The id 
electrons of iron experience a strong hybridization with 
the p states of oxygen and are widely spread over the 
whole valence band region. The total magnetic moment, 
calculated within the GGA, is 3.55 ps per f.u., which in¬ 
cludes the ferromagnetically induced O moments of about 
0.15 pB- 

The exchange interaction between Fe atoms is large 
and positive within the GGA: 20.5 meV and 2 meV be¬ 
tween the nearest and the next nearest neighbors, re¬ 
spectively. The estimated Curie temperature is about 
1000 K in the mean field approximation (700 K in the 
random phase approximation). It has to be mentioned 
here that the RPA provides usually a reasonable critical 
temperature, while the MFA represents average exchange 
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FIG. 2: (Color) The total and site-projected DOS of cu¬ 
bic BFO calculated within the GGA. The positive (negative) 
DOS values denote the majority (minority) spin states. 

interaction and overestimates Tq by 30 - 40%.^^ 

Thus, the Curie temperature, calculated in the cubic 
BFO using the GGA method, is several times larger then 
the experimental values (100 - 400 K). The reason is that 
the system is metallic within this approach and, there¬ 
fore, the double exchange interaction via the oxygens 
and the band ferromagnetism dominate over the other 
exchange mechanisms, which can reduce the magnetic 
interaction or change its sign. Therefore, the electronic 
correlations, which can be substantial in oxides, must be 
taken into account using a more appropriate method. 


2. GGA + U approach 

At the next step, we chose the GGA-I-U approximation, 
which provides a more adequate description of oxide ma¬ 
terials then the conventional LDA or GGA approaches.!^ 
Since we could not estimate the value of U from first prin¬ 
ciples, we used it as a parameter tracing the change of the 
electronic and magnetic properties with the value of U. 
In our simulations we applied the Hubbard U corrections 
on Fe 3d states and varied the value oi U^ff = U — J 
between 0 (here, GGA) and 9 eV. As expected, the occu- 
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FIG. 3: (Color) The total and site-projected DOS of cubic 
BFO calculated within the GGA-I-U method with (7e//=3 eV. 


pied (unoccupied) 3d Fe states are shifted down (up) in 
energy on the value of Ueff/2 that opens the gap in the 
minority spin channel for Ueff > 2 eV (see the DOS in 
Figj^with Ueff = 3 eV). The size of the spin band gap 
depends strongly on the Ueff value. 

In Fig|^ we show the total and all site-projected 
magnetic moments, which are plotted as a function of 
0 < Ueff < 9. The local magnetic moment of iron in¬ 
creases from 3 fiB to 4.18 fis within the given range 
of Ueff values, while the induced magnetic moment of 
oxygen decreases from 0.18 fiB to -0.08 fiB systemati¬ 
cally with the increase of Ueff- The magnetic moment 
induced on Ba is « 0.05 fj-B and its small value changes 
marginally on the model variations of U. For Ueff > 2 eV 
the total moment is integer (4 p-s/hu), which reflects the 
half-metallicity of the BFO. This is larger then the exper¬ 
imental values of 3.2 - 3.5 /rs/f.u., reported previously 
in Refs. I6I7I 

As next, we calculated the exchange interaction as a 
function of Ueff- The results are presented in Figj^ 
Since the BFO remains metallic within this approxima¬ 
tion, we can classify the main magnetic interaction as the 
double exchange. The consequence of the strong ferro¬ 
magnetic coupling is the largely overestimated Tc com¬ 
pared to the measured values. Although applying the 



FIG. 4: (Color) The total and site-projected magnetic mo¬ 
ments of cnbic BFO calculated as a function of U- Here, 
Ueff = 0 means the use of GGA. 


GGA-bU method reduces slightly the exchange interac¬ 
tion, mainly due to the reduction of the overlap between 
the Fe and oxygen states, the Gurie temperature is still 
too high in the whole range of Ueff- The dependence 
of Tc on Ueff is nontrivial. First, it reduces with the 
increase of Ueff, since the Fe 3d states are getting lo¬ 
calized and their overlap with sp states of oxygen. At 
Ueff = 2 eV, a band gap forms in the minority spin 
channel and the rise up with the increase of Ueff- 

Therewith, the exchange interaction between the nearest 
neighbors is getting larger only until Ueff « 4 eV and 
then decreasing, while the next nearest neighbor interac¬ 
tion continues to increase until Ueff ~ 6 eV. As result, 
the Gurie temperature gets up to 700 K within the RPA 
for Ueff « 6 eV and then starts to fall down. According 
to our estimations, the Gurie temperature agrees with ex¬ 
perimental values only for Ueff > 15 eV, which are non¬ 
physical for this compound. Thus, consideration only 
electronic correlation effects does not explain magnetic 
properties of BFO. This fact motivated us to improve 
the model of BFO and make further calculations. 


C. Simulation of XMCD spectra 

In complex transition metal oxides, the XAS and 
XMGD spectra at the L 2.3 absorption edges can be used 
as fingerprints of the ground state. In our study we sim¬ 
ulated the XAS and XMGD spectra varying the crystal 
structure, the chemical composition, the value of Ueff 
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FIG. 5: (Color) Exchange interactions (upper panel) and the 
Curie temperature (lower panel) of the cubic BFO calculated 
as a function of the Hubbard parameter Ueff applied on the 
Fe 3d orbitals. 


for electronic correlation effects, the magnetic order, and 
searching for optimal agreement with available experi¬ 
ments. 

Fig. [^displays the experimentally measurecP^ and the¬ 
oretically calculated Fe 2 ^ 2,3 XAS (upper panel) and 
XMCD (lower panel) spectra of the cubic BaFeOa. The 
La XAS spectrum exhibits a weak lower-energy shoulder 
at around 708 eV together with a high-energy shoulder 
at 713 eV. The L 2 XAS spectrum exhibits a lower-energy 
peak at around 721 eV. The Fe L 3 XMCD spectrum con¬ 
sists of major negative peak at 710 eV with low-energy 
shoulder and an additional positive peak at 713 eV. The 
Fe L 2 XMCD spectrum has double-peak structure. 

First, we found the best agreement with experiments 
was achieved for Ueff = 3 eV. This value of Ugff is 
used later for all our calculations if it is not mentioned 
specifically. However, the experimental results could not 
explained only by taking into account strong correlations 
effect. For the ideal cubic BFO structure, our simulations 
of the XAS and XMCD spectra provide x-ray absorption 
intensity in the Fe L 3 only at the major peak and high- 
energy shoulder, while the low energy peaks at the L 3 and 
L 2 edges is not reproduced (red dashed curve in Fig. |^. 



FIG. 6: (Color) Comparison between the experimentaP^ 
x-ray absorption (upper panel) and x-ray magnetic circu¬ 
lar dichroism (lower panel) at the Fe La.a edges of BaFeOs 
and the spectra calculated in the GGA-I-U approximation 
iUeff=3 eV) (black solid lines). Red dashed line show the 
spectra for ideal crystal structure, dotted blue lines present 
the theoretical spectra with oxygen vacancy. 


The calculation fails to reproduce low energy shoulder of 
the XMCD major negative peak at 709 eV in L 3 XMCD 
spectrum. Also, the simulations for the ideal cubic BFO 
phase provide only one a high energy peak structure in 
the L 2 XMCD spectrum, however, the experimental mea¬ 
surements exhibit a double-peak structure. 

The correct explanation of the experimental XAS and 
XMCD spectra is only possible by taking into account 
crystal imperfections, namely, oxygen vacancy. To inves¬ 
tigate the influence of oxygen deficiency on the XAS and 
XMCD spectra we created an oxygen vacancy in a dou¬ 
ble supercell in the first neighborhood of the second Fe 
atom. We found that the presence of the oxygen vacancy 
decreases the valency of nearest Fe ion from 4-|- to 3-1- 
(see the corresponding discussion below). The x-ray ab¬ 
sorption from the Fe atoms with the oxygen vacancy in 
close vicinity (dotted blue line in Fig. contributes to 
the low energy peaks in the L 2 absorption and XMCD 
spectra. The oxygen vacancy also explain the existence 
of the low-energy shoulder in Fe L 3 XMCD spectrum at 
709 eV. 

The same conclusion can be drawn from the interpre¬ 
tation of the XAS and XMCD spectra at the O K edge. 
Figj^ shows the experimentally measured x-ray absorp¬ 
tion and XMCD spectraP^ (open circles) at the O K edge 
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FIG. 7: (Color) Comparison between the experimental^ 
(circles) x-ray absorption (upper panel) and x-ray magnetic 
circular dichroism (lower panel) at the O K edges of BaFeOa 
and the spectra calculated by the LMTO method in the 
GGA-bU approximation ([/e//=3 eV) (black solid lines). Red 
dashed line show the spectra for ideal crystal structure, dotted 
blue lines present the theoretical spectra with oxygen vacancy. 


in BaFeOa together with theoretically calculated spec¬ 
tra for the ideal crystal structure of BaFeOa (red dashed 
curves) and the structure with oxygen vacancy (dotted 
blue curves). The calculations for the ideal structure are 
not able to describe the high-energy shoulder of major 
negative XMCD peak at 531 eV, which is due to the 
oxygen vacancy. We can conclude that the explanation 
of the experimental spectra at both the Fe La ,3 and O 
K edges demands the presence of crystal imperfections 
in BaFeOg. 


D. Fe valency and oxygen vacancies in BaFeOg-a 

To understand the vacancy formation we performed 
first-principles study of the Fe oxidation states in the 
cubic BaFeOg using a self-interaction method as it is de¬ 
scribed in Sec.|nj The Fe oxidation state vary in different 
compounds from 2+ to 4+. For a completely oxidized 
BFO, one can expect a high valency, Fe'^’*'. The total en¬ 
ergy study of the oxidation state and further comparison 
between BFO and BaFeOg-^ can provide us with some 
hints concerning the electronic properties and the prefer¬ 


able magnetic order in BFO. The Goodenough-Kanamori 
rules of the 180° superexchange suggest that the interac¬ 
tion Fe-O-Fe can be switched from antiferromagnetic to 
ferromagnetic when the Fe 3d-orbital filling changes. In 
this context, a variety of the superexchange options, such 
as Fe'^^-O-Fe"^’*', Fe^^-O-Fe^’*' and Fe'*+-0-Fe^+ need 
to be inspected. However, the Goodenough-Kanamori 
simplifications work robustly for insulators only. When 
theory deals with metallic magnetic oxide, such as the cu¬ 
bic BFO, the superexchange coupling does not dominate 
there and double exchange complicates the long-range 
interactions between magnetic species. 

In the case of the cubic symmetry, the Fe"^"*" electronic 
configuration is t^gCg, where the half-occupied Cg orbital 
has ether c? 322_^2 or representation. For the cu¬ 

bic BFO case, which was treated within the LSDA-SIC 
model, we found that the two Cg options result in the 
same equilibrium volume (see Fig^, which is substan¬ 
tially smaller then the experimental result (62.57 A^) and 
the volume calculated using the GGA (Fe^+). Although 
the volume computed within the GGA is very close to 
the experimental value, the corresponding total energy is 
about 22 meV higher then the total energies of Fe^"*" and 
Fe^'*' configurations. 

For the cation Fe^+, its t 2 g and eg levels must be fully 
occupied in the majority spin channel. This electronic 
configuration yields the total-energy minimum value, 
which is not far above that of the Fe^"*" ground state, 
as shown in Fig. The calculated equilibrium volume is 
in a better agreement with the experimental results then 
the values obtained for Fe^+ (SIG) and Fe^+ (GGA) con¬ 
figurations. The energetics suggest that both oxidations, 
Fe'^'*' and Fe^"*", can be realized in epitaxially fabricated 
thin films. Gonfigurations with other oxidation states are 
much higher in energy and are excluded from the further 
consideration. 

Further, we compared the total energies of Fe^'*' and 
Fe^"*" oxidation states in the presence of oxygen vacan¬ 
cies. For 6 = 0.03 the total energies in the equilibrium 
volumes are almost the same (see Fig. [^. Thus, the va¬ 
cancy formation lowers the energy of the Fe^"*" oxidation 
state, which agrees with our interpretation of the XAS 
and XMGD experiments. In a real BFO sample, both ox¬ 
idation states can coexist, which is as well evident from 
our XAS and XMGD simulations. In Fig. we present 
the total energy for a mixed configuration, which is lower 
then the others in a particular concentration range. 

As next, we investigated the impact of oxygen vacan¬ 
cies on magnetic properties of the BFO. The oxygen va¬ 
cancies were simulated within a 4 x 1 x I supercell, in 
which one oxygen site was occupied by an vacancy with 
a given probability. The disorder was simulated using 
the CPA method. The use of the supercell enables to 
take into account short-range effects in a simple man¬ 
ner, which can take place in real materials. Results for 
the exchange interaction, computed in this supercell and 
6 = 0.25, are presented in Fig. While the exchange in¬ 
teraction between the Fe"^"*" cations is strong and positive. 
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FIG. 8: (Color) The total energies of BaFeOa-^, calculated 
using the LSDA-SIC option for the Fe'*”'", Fe®"'", and a com¬ 
bination of them 1:1. Energies as functions of the volume 
without oxygen vacancies (a). Total energy as a function of 
volume and 5 (b). 



FIG. 9: (Color) Exchange interaction in the BaFe02.75 simu¬ 
lated in a 4 X 1 X 1 supercell. 
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EIG. 10: (Color) DOS of the BaEe02.75 simulated in a 4 x 1 x 1 
supercell. 


it is week and negative between the Fe^'*' atoms. The rea¬ 
son for this can be understood from the DOS, shown in 
Fig. Because of the vacancy (missing oxygen), the 3d 
states of the Fe^’*' are more localized (in particular the 
Cg states) and are shifted down in energy. This increases 
the tendency to an antiferromagnetic coupling. The cou¬ 
pling between the atoms in different oxidation states is as 
well antiferromagnetic and very strong due to the large 
difference of the DOS. This funding agrees with the main 
conclusions of S. Mori, who first developed a model for 
the magnetic structure of BaFeOa,.!^ 

Depending on the vacancy concentration and their 
distribution, the non-trivial exchange interaction in 
BaFeOa-a can results in various magnetic orders. In 
Fig. 11 we show the change of critical temperature with 
the oxygen vacancy concentration simulated within the 
MFA and RPA approaches. The negative Tq means here 
the absence of the ferromagnetic order (either antiferro¬ 
magnetic or non-collinear). Thus, for a moderate vacancy 
concentration S of 0.15 - 0.20, Tc lies in the experimental 
range of 250 - 100 K. 


IV. SUMMARY 

In this work we studied electronic and magnetic prop¬ 
erties of a cubic BaFeOa within a DFT framework using a 
GGA and a GGA-I-U approaches. We predicted that the 
BaFeOa in an ideal stoichiometric cubic structure should 
be half-metallic and ferromagnetic with a high Gurie tem¬ 
perature of 700 - 900 K. However, this finding disagrees 
the experiments, which found that the compound has 
a Tc far below the room temperature. To understand 
this discrepancy we simulated XAS and XMCD spectra 
comparing the theoretically obtained results with avail¬ 
able experimental data. Herewith, we varied the chemi¬ 
cal composition, structural and electronic parameters to 
obtain a better agreement with experiments. Our sim- 
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be explained by oxygen vacancies, which can be present 
in real BFO samples. Thereby, Fe atoms near an oxy¬ 
gen vacancy change their oxidation state from 4-1- to 3-1-. 
Thus, a real cubic BaFeOa may contain both Fe"^’*' and 
Fe^"*" cations, which explains the experimentally observed 
double peak structure in XAS and XMCD spectra. Fur¬ 
ther, we found that oxygen vacancies can substantially 
reduce the Curie temperature of a cubic BFO or even 
change the ferromagnetic order to an antiferromagnetic 
one or a non-collinear structure. Responsible for this is a 
strong antiferromagnetic exchange coupling between Fe 
atoms of different oxidation states. This result provides 
a clear route to fabricate a robust ferromagnetic BaFeOa 
with a high Curie temperature by reducing the amount of 
oxygen vacancies and stabilizing it in the cubic structure. 


FIG. 11: (Color) Critical temperature Tc in BaFeOa-j as a 
function of oxygen vacancy concentration 5. Here, the nega¬ 
tive Tc means the absence of the ferromagnetic order. 


ulations showed that the experimental observations can 
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